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We have investigated magnetic field dependences of a ferroelectric incommensurate-helimagnetic 
order in a trigonal magneto-electric (ME) multiferroic CuFei_ :E AL J 02 with x = 0.015, which ex- 
hibits the ferroelectric phase as a ground state, by means of neutron diffraction, magnetization and 
dielectric polarization measurements under magnetic fields applied along various directions. From 
the present results, we have established the H-T magnetic phase diagrams for the three principal 
directions of magnetic fields; (i) parallel to the c axis, (ii) parallel to the helical axis, and (iii) perpen- 
dicular to the c and the helical axes. While the previous dielectric polarization (P) measurements 
on 0uFei_ :E Ga x O2 with x — 0.035 have demonstrated that the magnetic field dependence of the 
'magnetic domain structure' results in distinct magnetic field responses of P [S. Seki et al, Phys. 
Rev. Lett., 103 237601 (2009)], the present study have revealed that the anisotropic magnetic field 
dependence of the ferroelectric helimagnetic order 'in each magnetic domain' can be also a source 
of a variety of magnetic field responses of P in CuFei_ x ^4 a: 02 systems (A =A1, Ga). 

PACS numbers: 75.80.+q, 75.25.+Z, 77.80.-e 



I. INTRODUCTION 

Nonlinear magneto-electric (ME) effects, in partic- 
ular magnetic-field control of ferroelectric polarization 
(P), in ME-multiferroics have been intensively studied 
since the discovery of the colossal ME-effcct in some 
magnetically frustrated transition metal oxides ,ii2i2di& 
In most of the ME-multifcrroics, the magnetic field de- 
pendences of P have been attributed to changes in the 
magnetic structures. For example, the magnetic-field- 
induced 90°-flop of P in TbMnC>3 has been explained 
by the first-order magnetic phase transition from the foe- 
plane cycloidal magnetic ordering to the ac-plane cy- 
cloidal magnetic ordering^ On the other hand, the 
recent experimental works have pointed out that mag- 
netic control of 'magnetic domain structures' can be 
an another ingredient for the ME-effects, in some rela- 
tively high-symmetry (trigonal, tetragonal or hexagonal) 
multiferroicsi^ i 10 ' 11 For example, Kimura et al. have 
argued that an anisotropic magnetic field dependence of 
P in a trigonal ME-multiferroic CuCr02 can be ascribed 
to the magnetic field dependence of the volume fractions 
of the three magnetic domains, in which the three mag- 
netic modulation wave vectors characterizing these do- 
mains are equivalent to each other because of the trigonal 
symmetry of the crystal structured This suggests that 
trigonal (or tetragonal, hexagonal) ME-multiferroics can 



provide opportunities for realizing a variety of magnetic 
field responses of P. 

In the past several years, a trigonal ME-multiferroic 
CuFe02 (CFO) has been the subject of increasing in- 
terest as a ME-multiferroic because of the discovery of 
the ferroelectricity in a magnetic-field-induced phased 
Subsequent studies have elucidated that the ferroelec- 
tric phase is stabilized even under zero magnetic field by 
substituting small amount of nonmagnetic Al 3+ or Ga 3+ 
ions for the magnetic Fe 3+ sitesi 13 ' 14 i 15 The magnetic 
structure in the ferroelectric phase has been determined 
to be a proper-screw-type helical magnetic structured 
The magnetic modulation wave vector is (q, q, |) where 
q = 0.202 ~ 0.210, and the helical axis is parallel to 
the [110] axis i 16 ' 17 Hereafter, we refer to the ferroelec- 
tric phase as ferroelectric incommensurate-magnetic (FE- 
ICM) phase. Recent polarized neutron diffraction studies 
have revealed that the spin-hclicity, left-handed (LH) or 
right-handed (RH) helical arrangement of spins, deter- 
mines the polarity of the local ferroelectric polarization 
emerging along the helical axis i 17 i 18 Because of the three- 
fold rotational symmetry about the hexagonal c axis (see 
FiglUa)), the magnetic ordering with the wave vector of 
(q, q, 1) results in three magnetic domains whose wave 
vectors of (q,q, §), (q, — 2q, |) and (— 2q, q, |) are crys- 
tallographically equivalent to each other, as illustrated 
in Fig. [He). I n this paper, we refer to these three do- 
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mains as 'g-domains'. 

Quite recently, Seki et al. have reported that the 'mag- 
netic domain structure' in the FE-ICM phase, specifi- 
cally the volume fractions of the three g-domains, can 
be controlled by applying a magnetic field in the trian- 
gular lattice planed They performed magnetization and 
dielectric polarization measurements on CuFei_ a; Ga 2 ;02 
(CFGO) with x = 0.035, in which the FE-ICM phase 
shows up as a ground state, under various magnitudes 
and directions of magnetic fields. These measurements 
revealed that the 'in-plane'-magnetic-field-dependence of 
the volume fractions of the g-domains results in distinct 
magnetic field responses of P, for example, 120°-flop of P 
by the magnetic fields rotating in the c-plane. However, 
magnetic field dependence of the magnetic ordering 'in 
each g-domain' was not discussed, because it is quite dif- 
ficult to extract the information on the magnetic ordering 
in each g-domain from the results of the macroscopic po- 
larization or magnetization measurements. In order to 
completely understand magnetic field dependences of P 
in slightly diluted CFO systems, therefore, it is critical to 
elucidate both of the magnetic field dependences of the 
'magnetic ordering in each g-domain' and the 'magnetic 
domain structure'. 

In the present study, we have investigated mag- 
netic field dependences of P and the FE-ICM order in 
CuFei_ x Al a; 02 (CFAO) with x = 0.015, which exhibits 
the FE-ICM phase as a ground state, by means of neu- 
tron diffraction, magnetization and dielectric polariza- 
tion measurements under various directions of magnetic 
fields. In order to elucidate the magnetic field depen- 
dence of the FE-ICM order 'in a g-domain', we have es- 
tablished H-T magnetic phase diagrams for the three 
principal directions of magnetic fields, specifically, (i) 
parallel to the c axis, (ii) parallel to the helical axis, and 
(iii) perpendicular to the c and the helical axes. The 
present results have revealed that the anisotropic mag- 
netic field dependence of the FE-ICM order 'in each q- 
domain' results in a variety of magnetic field responses of 
P as well as the magnetic field dependence of the 'mag- 
netic domain structure' does. 

We have also investigated a magnetic field dependence 
of sensitivity of P to a poling electric field (E p ). In 
our previous polarized neutron diffraction and in-situ 
pyroelectric measurements on CFAO(ir = 0.015) and 
CFGO (a; = 0.035) in zero field, we have found that the 
Al-substitution more significantly reduced the sensitiv- 
ity of P to E p as compared to the Ga-substitution. 17 
These measurements have also revealed that CFGO (a; = 
0.035) exhibits a homogeneous FE-ICM state, that is rel- 
atively close to the long-range-ordered state, in contrast 
to CFAO(a; = 0.015), that has an inhomogeneous domain 
state in zero fields We have thus concluded that the 'in- 
homogeneity' of the FE-ICM order, which must be rele- 
vant to the mobility of the magnetic domain walls, deter- 
mines the sensitivity of P to E p . On the other hand, our 
previous neutron diffraction measurements under applied 
magnetic fields showed that the homogeneous FE-ICM 




FIG. 1: (a) Crystal structure of CuFe0 2 . (b) The defini- 
tions of the hexagonal basis and the arrangements of the O 2- 
ions above (open blue circles) and below (filled blue circles) 
a Fe 3+ triangular lattice layer, (c) Schematic drawing of the 
q'-vectors of the three g-domains and the reciprocal lattice ba- 
sis, (d) Illustration of the magnetic structure in the FE-ICM 
phase and the directions of H« c , Hu q i and H± c q i. 



state can be realized even in CFAO(i£ = 0.015) by ap- 
plying a magnetic field along the c axisJ£ Therefore, the 
sensitivity of P to E p in CFAO(a; = 0.015) is expected to 
be controlled by an application of a magnetic field. This 
is the reason why we have selected CFAO(a; = 0.015) 
sample for the present study. 

This paper is organized as follows. In Sec. [TT1 we 
describe experimental details. Sec. IIIII consists of two 
subsections. In Sec. IIII Al we present anisotropic H-T 
magnetic phase diagrams for a g-domain using the results 
of the neutron diffraction measurements under the three 
directions of the magnetic fields, and also present the 
results of pyroelectric measurements under steady mag- 
netic fields. In Sec. IIII Bl we present magnetic field vari- 
ations of P and the FE-ICM order using the results of 
neutron diffraction, magnetization and dielectric polar- 
ization measurements under magnetic fields. In Sec. IIV1 
we summarize our conclusions. 
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II. EXPERIMENTAL DETAILS 

Single crystal of CuFei_ a: AL c 02 with x = 0.015 was 
prepared by the floating zone technique^. For the mea- 
surements of P, the sample was cut into a thin plate 
(~ 2 x 4 x 0.1 mm 3 ). Silver paste was applied on the 
[110] surface of the sample to form the electrodes. In the 
pyroelectric measurements under steady magnetic fields, 
the pyroelectric current was measured under zero electric 
field with increasing temperature, using an electrometer 
(Keithley 6517A). Before each pyroelectric measurement, 
we performed cooling with an applied poling electric field 
(E p ) from 15 K to 2 K. After the poling electric field was 
removed at 2 K, the sample was allowed to discharge for 
about 40 minutes in order to reduce residual current. The 
typical magnitude of E p in the present pyroelectric mea- 
surements was ~ 250 kV/m. We also investigated the 
E p dependence of P up to E p — 2.0 MV/m. External 
magnetic fields up to 5 T were provided by the Magnetic 
Property Measurement System (Quantum Design Inc.). 

We have performed dielectric polarization and mag- 
netization measurements on CFAO(a; = 0.015) under 
pulsed magnetic fields. The pulsed high magnetic fields 
up to 55 T were generated by a non-destructive magnet 
in the International MegaGauss Science Laboratory in 
ISSP, the University of Tokyo. The magnetization along 
the field direction was measured by the induction method 
using coaxial pick-up coils. The dimensions of the single 
crystal sample used for the magnetization measurements 
were 1.8 x 1.8 x 2.3 mm 3 . Following the pioneer work 
on the dielectric polarization measurement under pulsed 
magnetic fields by Mitamura et we detected the 

_ff-induced change in P by monitoring the polarization 
current through a voltage drop in the shunt resistance 
connected in series. By integrating the current with re- 
spect to time, we obtained the magnetic field variations 
of P. In this measurement, the poling electric fields (E p ) 
were continuously applied during all the measurements 
of the P-H curve. The magnitude of E p was typically 
^200 kV/m. The single crystal sample for this measure- 
ment was identical to the sample used for the present 
pyroelectric measurements. 

The neutron diffraction measurements under applied 
field were carried out at the two-axis neutron diffrac- 
tometer E4 installed at Berlin Neutron Scattering Cen- 
ter (BENSC) in Helmholtz Centre Berlin for Materials 
and Energy. A typical dimension of the single crystal 
samples for the neutron experiments was ~3x4x 4.5 
mm 3 . External magnetic fields along the hexagonal [001], 
[110], and [110] directions were provided by the cryomag- 
nets HM-1, HM-2 and VM-1, whose maximum fields are 
6T, 4T, and 14.5T, respectively. Note that the mag- 
netic field along [001] axis was canted by ~ 12° from 
the [001] axis toward [110] axis, because of limitations 
of the windows of the horizontal field cryomagnet HM-1. 
For all the neutron diffraction measurements, the sample 
was mounted in the cryomagnet with a (h, h, I) scattering 
plane. In the experiments with the [001] and [110] mag- 



netic fields, the collimation was 40 '-40 '-40', and a sin- 
gle detector was used. In the experiment with the [110] 
magnetic field, the collimation was 40'-40'-open, and a 
2D-position-sensitive-detector was used. 

In this paper, we have employed a conventional hexag- 
onal basis defined as shown in Fig. [TJa), in order to 
describe the directions of the magnetic propagation vec- 
tors of the three g-domains, although structural transi- 
tions from the original trigonal structure to a monoclinic 
structure have been reported for some of the magnet- 
ically ordered phases (including the FE-ICM phase) of 
CuFeO^^iZi and CuFei-^AUO^. 



III. RESULTS AND DISCUSSIONS 

A. Magnetic field dependence of the phase 
transitions in a g-domain 

1. Neutron diffraction measurements under steady 
magnetic fields 

We firstly investigated the temperature variations of 
the magnetic ordering in a g-domain under steady mag- 
netic fields, by means of the neutron diffraction measure- 
ments. In order to define the relationship between the di- 
rections of the magnetic fields and the magnetic structure 
in a g-domain, we introduce the c-plane-projection of the 
q- vector, g' = (g,g, 0). In the present neutron diffraction 
measurements, we applied magnetic fields along the three 
directions; (i) nearly parallel to the [001] axis,- 5 - (ii) par- 
allel to the [110] axis, and (iii) parallel to the [110] axis. 
Since the magnetic reflections on the (h, h, I) scattering 
plane belong to the g-domains with the wave vector of 
(<7)9, §), the directions of the [110] axis corresponds to 
the g'-vector of the (g, g, |)-domain. Hereafter, we refer 
to the magnetic fields along these three directions as H\\ c > 
H\\ q < and H ±c y (see Fig. QJd)). 

Before discussing the present results, we should review 
the magnetic phase transitions in CFAO (x = 0.015) 
in zero magnetic field. As reported by Terada et al.^ 
CFAO(x = 0.015) has three magnetically ordered phases 
in zero magnetic field. The typical magnetic diffraction 
profiles of the (h, h, |) reciprocal lattice scans in zero field 
are shown in Fig. H]by black filled circles. With decreas- 
ing temperature from the paramagnetic (PM) phase, 
two collinear-incommensurate magnetic phases show up; 
the higher temperature phase is the oblique-partially- 
disorderd (OPD) phase and the lower temperature phase 
is the partially disordered (PD) phase. Both of the mag- 
netic reflections corresponding to the OPD and PD mag- 
netic orderings are assigned as (g, g, |). The incommen- 
surate wave number for the OPD phase, gopD ~ 0.195 
is almost independent of temperature, while that for the 
PD phase, gpo, varies from 0.20 to 0.22 with decreas- 
ing temperature. Previous neutron diffraction measure- 
ments by Terada et a/. 26 i 27 have revealed that the mag- 
netic moments in the PD and OPD phases are canted 
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FIG. 2: [(a-l),(b-l),(c-l)] The relationships between the applied magnetic fields and g'-vector of the g-domain observed in 
the present neutron diffraction measurements under (a-1) Hu c , (b-1) Hu q i and (c-1) H ±cq i. The magnetic diffraction profiles 
measured in the cooling processes under [(a-2)-(a-4)] Hu c , [(b-2)-(b-4)] H\\ q , and [(c-2)-(c-4)] H± cq /. The data observed in the 
magnetic fields are shown by the open circles, and the those in zero field are shown by the black filled circles. Note that the 
difference in the width of the magnetic diffraction profiles between the measurements under H« q i and the others is due to the 
different experimental configurations, specifically the collimations and the detectors. See Sec(IT]for details. 



by about 12°, 50° from the [001] direction toward the 
[110] direction, respectively. The ground state of CFAO 
(x = 0.015) is the FE-ICM phase, as mentioned in in- 
troduction. The magnetic diffraction profile in the FE- 
ICM phase is characterized by the two magnetic reflec- 
tions assigned as (q,q, f ) and (| — q, | — q, using 
the hexagonal basis, as shown in Fig. [7[a-3). The mag- 
netic modulation wave number in the FE-ICM phase, 
9FE-ICM = 0.202 ~ 0.210, slightly depends on the Al- 
concentration and applied magnetic fields, but is almost 
independent of temperature. It should be noted that 
a small peak at (|,-j,§) corresponds to the collinear- 
commensurate 4-sublattice phase, which coexists with 
the FE-ICM phase at low temperatures because of a 
slight macroscopic inhomogeneity of the Al-concentration 
in the samples— & 

We now discuss the results of the present neutron 
diffraction measurements under applied magnetic fields. 
In Fig. [21 we show the magnetic diffraction profiles 
measured on cooling under H» c , H^ ql and H± c ^ q i of 
4 T. In all the three cooling, the successive magnetic 
transitions [OPD— >PD^FE-ICM] were observed, and 
no strong magnetic field dependences of the magnetic 
diffraction profiles were found in the OPD and PD phase, 



although the wave number of the PD phase is slightly de- 
pendent on -ff|| c . However, we found that the magnetic 
diffraction profiles in the FE-ICM phase under H\\ c , H\\ q i 
and H± c , q > are remarkably different from each other, as 
shown in Figs. IKa-3), EHb-3) and[i;c-3). 

Figure [U[b-3) shows the magnetic diffraction profile in 
the FE-ICM phase under H\\ q t, suggesting that the mag- 
netic field applied along the helical axis does not result 
in the drastic change in the FE-ICM order. In contrast, 
the magnetic fields applied perpendicular to the helical 
axis, _ff|| c and H± c ^ q i, significantly affect the magnetic 
diffraction profiles in the FE-ICM phase. Figure [2Ja-3) 
shows that H\\ c sharpens the magnetic diffraction profile 
in the FE-ICM phase, as was partly reported in Reffl6l. 
This suggests that the FE-ICM order under H\\ c is rel- 
atively close to the long-range-ordered state, while the 
FE-ICM order is inhomogeneous domain state in zero 
magnetic field^ as mentioned in introduction. This 
_ff|| c -dependence of the magnetic correlation might be at- 
tributed to the local lattice distortions due to the Al- 
substitution. 

The previous studies on CFAG±£ and CFGG±£i±£ have 
pointed out that the Al-substitution should result in lo- 
cal lattice distortions because the ionic radius of Al 3+ 
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FIG. 3: The temperature variation of the magnetic diffraction 
profile of the (h, h, |) reciprocal lattice scans under H ±cq i 
of 14.5 T. Inset shows relationship between the directions of 
H± Ct gi and the magnetic moments of the PD magnetic order 
in zero field. 



is much smaller than that of Fe 3+ , and have also sug- 
gested that this local lattice distortion strongly disturbs 
a coherent magnetic ordering in the FE-ICM phase be- 
cause the distortions randomly lift the local degeneracy 
in the competing exchange interactions. On the other 
hand, the previous synchrotron radiation x-ray diffrac- 
tion measurements on pure CFO under applied magnetic 
fields revealed that the lattice constants in the FE-ICM 
phase vary with Hn c \ 24 ' specifically, the lattice constant 
along the [1 10] direction, which is the b axis in the mon- 
oclinic notation, linearly decreases with increasing Hu c . 
Taking account of these results, one can expect that the 
local lattice distortion relaxes with increasing Hu c . For 
further investigation on this problem, x-ray diffraction 
measurements on CFAO system under applied magnetic 
fields are required. 

In contrast to Hn c , in the cooling process under H± cq i, 
the magnetic diffraction corresponding to the FE-ICM 
order is rather diffusive. This implies that the FE-ICM 
order is suppressed by -Hj_ Cj9 '. Actually, as shown in Fig. 
|31 the PD to FE-ICM magnetic phase transition is not de- 
tected in the cooling process under H± q i, c = 14.5 T, and 
instead, the PD magnetic ordering survives even at 2.9K. 
This suggests that -Hj_ c , 9 ' favors the PD magnetic order- 
ing, whose magnetic moments lie nearly perpendicular 
to the magnetic field, rather than the proper-screw-type 
magnetic order in the FE-ICM phase. 



2. H-T magnetic phase diagram for H« 



H\\ q r and H± cq i 



In Figs. H{a)-Uc), we now present the H-T magnetic 
phase diagrams for H\\ c , H\\ q i and H± Ctq > deduced from 
the present results. Although the present neutron diffrac- 
tion measurements could not reach the high field region 
of the H\\ q i-T phase diagram, we have confirmed that 
the FE-lCM phase extends to the high field region by 
the present dielectric polarization measurements under 
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FIG. 4: [(a)-(c)] The H-T magnetic phase diagrams for (a) 
ffiic! (b) H\\ q i and (c) H± c q i. Open circles and solid lines de- 
note the phase boundaries determined by the present neutron 
diffraction measurements on cooling. The phase boundaries 
in the high field region of the H« c -T phase diagram (dashed 
lines) were drawn in analogy of the Hu c ~T phase diagram of 
CFAO(a; = 0.02)>2£ To check the reasonability of the phase 
boundaries in (a), we have shown the transition field between 
the FE-ICM phase and the FI phase determined by the dielec- 
tric polarization measurements under pulsed magnetic fields 
presented in Sec. HUB II (open triangles). 



pulsed magnetic fields up to 30 T to be mentioned in Sec. 
IIIIB II Here, we should emphasize that these phase di- 
agrams represent the magnetic ordering 'in a g-domain', 
and moreover, those revealed that the FE-ICM order in 
a g-domain shows the anisotropic responses for the in- 
plane magnetic fields of Hu q i and H± c _ q r, which was not 
directly observed in the previous macroscopic polariza- 
tion and magnetization measurements p^i 



3. Pyroelectric measurements in steady magnetic fields 

We have also performed pyroelectric measurements 
under applied magnetic fields, to observe the magnetic 
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FIG. 5: [(a-l),(b-l),(c-l)] The relationships between the applied magnetic fields and g'-vectors for (a-1) H\\ c , (b-1) H\\e~ and 
(c-1) H± Ct E p - The magnetic field dependence of the temperature variations of P under (a-2) H\\ c , (b-2) H\\ Ep and (c-2) Hj_ C: e p - 
(d) The definition of the H±/ E for the in-plane magnetic fields of H^ Ep and Ha_ c ,e p - (e) The H^J E -dependences of P at 2.0 
K normalized to the value of P in zero field. The dashed line is the guide to eyes. 



field dependences of P corresponding to the FE-ICM or- 
der observed in the present neutron diffraction measure- 
ments. As mentioned in introduction, the direction of 
the spontaneous electric polarization in the FE-ICM or- 
der is parallel to the helical axis,— that is, the direc- 
tion of the g'-vector. In addition, the previous polarized 
neutron diffraction measurement under applied electric 
field revealed that an application of a poling electric field 
cannot result in a single g-domain state . 17 ' 18 Therefore, 
the measured electric polarization has to be the sum of 
the contributions from the domains with the different q- 
vectors. This situation prevent us from investigating the 
anisotropic ME-responses in a q-domain by pyroelectric 
measurements . 

To overcome this problem, we have thus selected the 
[120] plane, which is crystallographically equivalent to 
the [110] and [210] planes, as the electrode surfaces, and 
applied magnetic fields along three principal directions; 
(i) parallel to the c axis, (ii) parallel to the poling electric 
field (E p ), and (iii) perpendicular to the c axis and E p . 
We refer to these three directions of the magnetic fields 
as H\\ c , H\\ Ep and H± CjEp , respectively. The directions 
of the poling electric field, the g'-vectors of the three do- 
mains and the applied magnetic fields are schematically 
drawn in Figs. [3 a-1), ODJb-l) and[5^c-l). In this config- 
uration of the electrodes, only the two g-domains with 
the wave vectors of (g, g, |) and (g, — 2g, |) contribute to 
the measured electric polarization, because the electric 



polarization vector in the (— 2g,g, |)-domain is perpen- 
dicular to the normal vector of the electrode surfaces. In 
addition, the H\\e p ~ (or H± c ,e p -) dependence of the FE- 
ICM order in the (g,g, |)-domains is expected to be the 
same as that in the (g, — 2g, |)-domains, because of the 
symmetry of the crystal and magnetic structures. By this 
configuration of the electrodes, the anisotropic in-plane- 
field-dependences of P can be observed. 

Figures GHb-2) and GDJc-2) show the temperature vari- 
ations of P under H\\e p an d H± c .e p - For both of -ff||_E p 
and H± C ,E , the magnitude of P decreased with increas- 
ing magnetic field. Taking account the present neutron 
diffraction measurements revealing that the FE-ICM or- 
der was significantly affected by H± c _ q > and was less af- 
fected by -ff|| g ', we conclude that the H± Ctq r -components 
of the magnetic fields are relevant to the reduction of P. 
Actually, H± c ^e p more remarkably reduced P than H\\e p - 
We thus show, in Fig. EJe), the values of P at T = 2 K 
normalized to the values in zero magnetic field as a func- 
tion of the effective magnetic field applied perpendicular 
to the c axis and the g'-vector, H e /J q , (see Fig. EJd)), 
specifically 

H ±/, q > = H ±c,E p cos 30° = H\\ Ep sin 30°. (1) 

We found that the H± Ct E p - and -ff||£ p -dependences of P 

at T = 2.0 K are scaled by H^/ q/ , This clearly shows 
that Hj_ Ctq i dominates the 'in-plane' magnetic field de- 
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(a-1) T=6.0K,H|| [110] 



(b-1) T= 2.8K, H|| [110] 
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FIG. 6: The ^-dependences of the values of P at T = 2.0 
K in zero field, under H\\ c = 5.0 T and H ±c>Ep = 5.0 T. The 
inset shows the magnification of the region of < E p < 300 
kV/m. The solid gray lines are guides to eyes. 



pendence of P, as was expected from the results of the 
present neutron diffraction measurements. 

As for the origin of the reduction of P under a mag- 
netic field having the H± Ciq i -component, we can propose 
three possibilities; the first one is that the volume frac- 
tion of the FE-ICM order in each g-domains was reduced 
because of the PD magnetic order retained by the mag- 
netic field, the second is that the magnetic structure of 
the FE-ICM order was partly modified by the magnetic 
field, and the third is that the sensitivity of P to E p 
was reduced because of the reduction of the coherence 
of the magnetic order in the FE-ICM phase under the 
magnetic field. A combination of two or three of them 
is also possible. Although we cannot identify the ori- 
gin of the reduction of P only from the present results, 
but the i?p-dependence of P in H± c ^e p = 5 T shown in 
Fig. [6] suggests that the third scenario contributes to the 
reduction of P. 

In Fig. [5Ja-l), we show the temperature variations 
of P under Hu c . In contrast to the in-plane magnetic 
fields, the application of Hn c enhances P. Moreover, 
Fig. [6] shows that the sensitivity of P to E p was sig- 
nificantly enhanced by applying H» c . As mentioned in 
introduction, the previous study on CFAO and CFGO 
has revealed that the sensitivity of P to E p is deter- 
mined by the 'inhomogeneity' of the FE-ICM order, that 
must be relevant to the mobility of the magnetic do- 
main walls in the FE-ICM phase. 17 Taking account of 
the present neutron diffraction measurements revealing 
that CFAO(a; = 0.015) exhibits the long-range-ordered 
FE-ICM phase above Hu c — 4 T, we attributed to the 
enhanced sensitivity of P to E p to the i?|| c -dependence 
of the magnetic correlation in the FE-ICM phase. 
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FIG. 7: [(a),(b)] The //x c ,ij'-variations of the neutron diffrac- 
tion profiles of the (h, h, |) reciprocal lattice scans at (a) 6.0 
K and (b) 2.8 K. [(a-2), (a-3), (b-2), (b-3)] Schematic drawings 
of the H± ,g' -induced magnetic phase transitions in each q- 
domain. The sizes of the arrows qualitatively show the volume 
fractions of the g-domains. 



Note that the curvatures of the E p dependence of P 
implies that the saturation value of P was also enhanced 
by H\\ c - This might be ascribed to the il|| c -dependence of 
the (local) magnetic structure including the wave number 
of the magnetic order in the FE-ICM phase. 



4- Neutron diffraction measurements under applied 
magnetic field 

On the basis of the results of the field-cooling scans, 
we now discuss magnetic field variations of the FE-ICM 
order including those of the domain structure. 

According to the H± c ^ q i-T phase diagram shown in Fig. 
2{c), the magnetic phase transition from the noncollincar 
FE-ICM phase to the collinear PD phase is expected in 
a -ffj_c, (/-increasing process. Figures EI a-1) and [^b-l) 
show the iJj_ Ciq '-variations of the neutron diffraction pro- 
files measured at T = 6.0 K and 2.8 K, respectively. In 
the magnetic field scan at T = 6.0 K, the intensities 
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T= 6.0 K, H|| C = 0T-^6T-^ 0T 

g 




0.19 0.21 0.23 0.27 0.29 0.31 
(h,h,3/2) 



FIG. 8: The ff|| c -variation of the neutron diffraction profiles 
of the (h, h, |) reciprocal lattice scans at T = 6.0 K. 



of the magnetic reflections corresponding to the FE-ICM 
order monotonically decrease with increasing H± c q i, and 
disappear around 9.0 T, as expected from the phase dia- 
gram. On the other hand, magnetic reflections described 
by the wave vector of the PD phase, (<?pd,9pd, §) with 
<Zpd ~ 0.22, emerge above 8.0 T. This result apparently 
manifests the magnetic field induced phase transition 
from the FE-ICM phase to the PD phase. 



B. Magnetic field variations of the FE-ICM order 

In the magnetic field scan at T = 2.8 K, the field in- 
duced magnetic transition was also observed, as shown 
in Fig. Efb-1). However, we found that the diffraction 
profile of the field-induced PD phase is rather diffusive. 
This indicates that the sinusoidally amplitude modulated 
magnetic structure of the PD phase is no longer stable 
at low temperatures. In addition, we also found that the 
intensity corresponding to the field induced PD magnetic 
order is considerably small. This implies that the volume 
fraction of the (q, q, |)-domains is reduced by the applied 
magnetic field, namely that the volume fractions of the 
three g-domains are changed by -Hj_c,<j', as illustrated in 
Figs. [DJb-2) and[7]Jb-3). This indicates that at low tem- 
perature, the in-plane magnetic field favors the proper 
helical magnetic domains whose helical axis follows it, as 
demonstrated in the previous worksi 11 ' 23 



To summarize the -ffj_ Cjl3 '-dependence of the magnetic 
ordering in this system, the magnetic phase transition 
from the FE-ICM phase to the PD phase is observed at 
the relatively high temperature (T = 6.0 K), and the 
re-population of the g-domains occurs at low tempera- 
tures (T = 2.8 K). As mentioned in introduction, Seki et 
al. have demonstrated that a magnetic field rotating in 
the triangular lattice plane can induce 120°-flop of the 
electric polarization because of the re-population of the 
g-domainsj^i The maximum magnitude of the rotating 
magnetic field shown in Ref. 11 was 6.5 T, and the mea- 
surements were carried out at 2.0 K. However, the present 
results suggest that more diversity in the ME-responses 
should be found by applying magnetic fields beyond ~ 9 
T and by changing the temperature. 

Figure[5]shows the i?|| c -variation of the neutron diffrac- 
tion profiles at T = 6.0 K. As seen in the field-cooling 
scans, the magnetic field applied along the c axis sharp- 
ens the magnetic diffraction profile. In addition, the 
wave number of the FE-ICM order, which was distributed 
around q ~ 0.21 in zero field, was concentrated at 
q = 0.207 above Hn c = 4T. We also found that the sharp 
magnetic diffraction profiles retained after removing the 
magnetic field. This history dependent behavior indi- 
cates the long-range-magnetic ordering realized by the 
application of H\\ c remained even after returning to the 
zero field. 



1. magnetization and dielectric polarization measurements 
under pulsed magnetic fields 

We also performed the dielectric polarization and mag- 
netization measurements under pulsed magnetic field up 
to 30 T. The electrode configuration and the magnetic 
field directions were selected to be the same as those in 
the pyroelectric measurements discussed in Sec. 1111 A 31 

We firstly discuss the metamagnetic transition in this 
system using the results of the magnetization measure- 
ments. In H\\ c up to 30 T, three magnetic phases ap- 
pear, as shown in Figs. [9ja-2)-[9ja-3). In the previous 
work on CFAO(a; = 0.02) ^ Terada et al. have reported 
that the first field induced phase from the FE-ICM phase 
is the slightly incommensurate magnetic phase referred 
to as the 'FI phase'. We thus considered that the first 
field induced phase from the FE-ICM phase in CFAO 
{x = 0.015) is the FI phase. Around if|| c = 20 T, where 
the induced magnetization approaches ~ 5/3/xb, the sys- 
tem undergoes further magnetic phase transition. In the 
in-plane magnetic fields, magnetization plateaus with the 
magnetization of ~ 5/3/ib were also observed around 
23 T. Comparing the results of the previous magneti- 
zation measurements on pure CFO under pulsed mag- 
netic fields we realize that this magnetic phase is the 
3-sublattice (3SL) phase, which is also observed in pure 
CFO^ 

Let us move on to the results of the dielectric polar- 
ization measurements. In H\\ C1 the finite electric po- 
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FIG. 9: [(a-l),(b-l),(c-l)] The relationships between the directions of E p , the g'-vectors and the applied magnetic fields of (a-I) 
i/|| c , (b-I) H\\e p and (c-f) H± Ct E p - The results of the polarization and magnetization measurements under the pulsed magnetic 
fields of [(a-2),(a-3)] H\\ c , [(b-2),(b-3)] H\\ Ep and [(c-2),(c-3)] H± CyEp . [(d)-(e)] The i/|| c -dependence of the magnetization and 
P in CFAO(a; = 0.015) up to 55 T. The inset shows the ifii (.-dependence of dMjdH\\ c around the fifth-field-induced phase 
transition. 



larization was observed only in the FE-ICM phase, as 
shown in Figs. [9^a-2)-[9^a-3). In the 7J|| c -increasing pro- 
cess, P rapidly increases in the magnetic field region 
of 5 < i?|| c < 12 T. This enhancement corresponds to 
the magnetic field dependences of the magnetic correla- 
tion and the magnetic modulation wave number, which 
were observed in the present neutron diffraction mea- 
surements. In the 7J|| c -decreasing process, P is also ob- 
served to emerge in the FE-ICM phase. In contrast to the 
-ff|| c -increasing process, the value of P is almost indepen- 
dent of the magnetic field, and is larger than the value 
in the zero field state before the measurement. Judging 
from the symmetry of the magnetic and crystal struc- 
ture in this system, it is reasonable to assume that the 
volume fractions of the three g-domains are not changed 
by .ff|| c . Hence, we ascribed this iJ|| c -variations of P to 
the history dependent behavior of the FE-ICM order in 
a g-domain observed in the present neutron diffraction 
measurements. 

Figures [9jb-3) and[9]Jc-3) show the 'in-plane' magnetic 
field dependences of P measured at relatively low tem- 
perature, T — 2.0 K. From the results of the present 
neutron diffraction measurements, it is expected that the 
magnetic-field-induced re-populations of the q-domains 
result in anisotropic magnetic field variations of P. How- 
ever, the expected anisotropic behaviors are not ob- 
served, although the details of the P-H± Ct E p and P-H\\e p 



curves are slightly different from each other. This result 
suggests that the sweeping rate of the pulsed magnetic 
fields were too fast to induce the re-population of the 
q- domains. 

On the other hand, at T = 6.0 K, we found that 
the in-plane field variations of P were quite anisotropic, 
as shown in Figs. E)^b-2) and E^c-2). In the H± c ^e p - 
increasing process, the value of P started to decrease 
around 8 T, and disappeared around 15 T, which is far 
below the transition field to the 3SL phase. In contrast, 
in the H\\e p -increasing process, the finite value of P is ob- 
served up to the transition field to the 3SL phase. These 
results revealed that the 'in-plane' field induced FE-ICM 
to PD transition in a q-domain occurs even in the pulsed 
magnetic fields. It should be noted that this field induced 
FE-ICM to PD phase transition is not clearly observed 
in the magnetization measurement. This is because the 
FE-ICM to PD phase transition occurs only in the q- 
domains with the wave vectors of (q, q, |) and (— q, 2q, |) 
while all the three q-domains contributes to the measured 
magnetization. 

We also found that in the H± c ^e p - and iJ|| Bp -decreasing 
process, the values of P are significantly smaller com- 
pared to the values in zero field before the measurements, 
in contrast to the -ff|| c -decreasing process, in which P was 
comparable to the value at the initial state. At this stage, 
we have no clear explanation for this result. In order to 
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investigate this point, further neutron diffraction mea- 
surements under applied magnetic fields are required. 

In Figs. M,d) and[9|e), we show the H|| c -dependence 
of the magnetization and P up to 55 T, respectively. 
Although, in this paper, we do not focus on the high- 
magnetic-field behavior in this system, it is worth 
mentioning here that the fifth-magnetic-field-induced 
phase transition, which has been recently found in pure 
CFOr^£ was also observed in CFAO(x = 0.015) around 
H = 48 T. No finite electric polarization was detected 
in the magnetic-field-induced phases including the 'fifth' 
phase. 

IV. CONCLUSION 

We have investigated the magnetic field dependence 
of the FE-ICM order in the trigonal ME-multiferroic 
CuFei-zALrC^ with x = 0.015 by means of the neu- 
tron diffraction, dielectric polarization and magnetiza- 
tion measurements under the magnetic fields applied 
along various directions. 

We have established the H-T magnetic phase diagrams 
for the three principal directions of the magnetic fields, 
Hu c , H\\ q i and H± Ctq >. It should be emphasized that 
these H-T phase diagrams represent the magnetic or- 
dering 'in a g-domain', and reveal the anisotropic in- 
plane-magnetic-field responses of the FE-ICM order in 
a q-domain, which was not directly observed in the pre- 
vious macroscopic (bulk) polarization and magnetization 
measurements^ 

We have also found that the sensitivity of P to E p 
in this system was controlled by applying a magnetic 



field. This indicates that the 'inhomogeneity' of the 
FE-ICM order, which is also controlled by an applied 
magnetic field, is relevant to the sensitivity of P to 
E p . This result is consistent with our previous polar- 
ized neutron diffraction study on CFAO(a; = 0.015) and 
CFGO(x = 0.035)^ 

While the recent dielectric polarization measurements 
on CuFei_a;Gaa;02 with x = 0.035 by Seki et al. have 
demonstrated that the magnetic field dependence of the 
'magnetic domain structure', specifically the volume frac- 
tions of the three g-domains, results in the distinct mag- 
netic field responses of Pr^ the present results have re- 
vealed the anisotropic magnetic field dependence of the 
FE-ICM order 'in each q-domain' can be also a source of 
a variety of the magnetic field dependence of P in this 
system. 
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